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ABSTRACT 

We consider a model in which the ultra-relativistic jet in a gamma-ray burst (GRB) is cold 
and magnetically accelerated. We assume that the energy flux in the outflowing material is 
partially thermalized via internal shocks or a reverse shock, and we estimate the maximum 
amount of radiation that could be produced in such magnetized shocks. We compare this es- 
timate with the available observational data on prompt y-ray emission in GRBs. We find that, 
even with highly optimistic assumptions, the magnetized jet model is radiatively too ineflicient 
to be consistent with observations. One way out is to assume that much of the magnetic en- 
ergy in the post-shock, or even pre-shock, jet material is converted to particle thermal energy 
by some unspecified process, and then radiated. This can increase the radiative efficiency suf- 
ficiently to fit observations. Alternatively, jet acceleration may be driven by thermal pressure 
rather than magnetic fields. In this case, which corresponds to the traditional fireball model, 
sufficient prompt GRB emission could be produced either from shocks at a large radius or 
from the jet photosphere closer to the center 

Key words: acceleration of particles - MHD - radiation mechanisms; non-thermal - rela- 
tivistic processes - shock waves - gamma-ray burst: general 



1 INTRODUCTION 

A great deal of progress has been made in our understanding 
of gamma-ray bursts (GRBs), thanks to the launch of a number 
of dedicated satellites (BeppoSAX, HETE-2, Swift and Integral). 
These satellites rapidly communicate burst locations to ground- 
based optical and radio telescopes, which has enabled detailed fol- 
low up study of the GRB afterglow emission. It is now known 
that GRBs produce highly relativistic and beamed jets contain- 
ing energy ~ 10^' erg (s ee Meszaros' '2002]; IPiranI 120051: IZhand 
1 20071 : iGehrels. RamirezT^ uiz & Fox 200f^ for extensive reviews 
of these and other developments). It is also well established that 
there are two classes of GRBs. One class, called long-GRBs — 
those lasting for more than a few seconds — is produced when a 
massive star collapses a t the end of its nuclear burning life (see 
IWooslev & Bloomll200^ for a review). For the other class, called 
short-GRBs - those lasting for less than a few seconds - at least 
some members are believed to result from mergers of compact stars 
in binary systems (Gehrels et al. 2009, and references therein). 

Despite this impressive progress, several fundamental ques- 
tions remain unanswered. Foremost among these is the composition 
of the relativistic jets that power GRBs. We do not know whether 
GRB jets consist of a normal proton-electron plasma or if they are 
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dominated by electron-positron pairs. Furthermore, it is uncertain 
whether the jets are dominated by matter or magnetic fields (Poynt- 
ing outflow). The related question of how the observed y-ray radi- 
ation is produced is also poorly understood. 

A popular model for converting jet energy to particle ther- 
mal energy and thereby to radiation is the internal shock model 
(Naravanetal. 1992; Rees, Meszaros 1994; Sari & Piran 1997*). 
According to this model, the relativistic wind from the central en- 
gine of a GRB has a variable Lorentz factor, which leads to colli- 
sions between faster and slower moving ejecta. A fraction of the ki- 
netic energy of the jet is converted to thermal energy in these "inter- 
nal" shocks. A fraction of this thermal energy then goes into elec- 
trons and is rapidly radiated away as y-ray photons via synchrotron 
and inverse-Compton processes. The internal shock model natu- 
rally produ ces the rapid vari ability observed in the y-ray emission 
of GRBs (Sari fePiranll 1997 1. This is one of its principal virtues. 

The internal shock model, however, has a problem, viz., 
the efficiency (see eq. [12] for the definition ) for convert- 
ing jet energy to radiation is relatively low (■Kumar '19991 
. Lazzat i. Ghisellini & Celotti 1999; Panaitescu, Spada & Meszar(3 
1 19991) . The efficiency depends on the relative Lorentz factor of the 
colliding blobs, and also, in the case of magnetized ejecta, on the jet 
magnetization parameter cr (defined in eq.[2j. Since the efficiency 
6y of a GRB can be measured directly from observations of the 
prompt and afterglow emission, we can constrain the parameters 
of the internal shock model, notably the magnetization cr of the jet 
material. 
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Another location where the jet energy may possibly be con- 
verted to y-rays might be the deceleration radius where the jet starts 
to slow down as a result of its interaction with an external medium. 
Two shocks are formed in this interaction, one of which, the "for- 
ward" shock, heats up the external medium and produces the af- 
terglow emission, and the other, the "reverse" shock, propagates 
into the GRB jet. The energy released in the reverse shock could 
be radiated as y-rays via synchrotron and/or inverse-Compton pro- 
cessesQ The efficiency for converting jet energy to y-rays depends 
on various parameters, including the cr of the jet material. 

Thus, in either the internal shock model or the reverse 
shock model, the y-ray efficiency 6^ depends on the magne- 
tization (T of the jet. The present work is motivated by the 
fact that, under some circumstances, we can independently es- 
ti mate cr for a GRB jet. This follow s from the recent work 
of iTchekhovskov. Naravan & McKinnevI j2010b ) who studied the 
properties of a magnetically accelerated GRB jet. If the jet material 
is cold, i.e., there is no thermal pressure, and all the acceleration 
is from electromagnetic forces (Poynting-dominated jet), these au- 
thors show that cr can be estimated from the terminal Lorentz factor 
yj and the opening angle dj of the jet. Both of the latter quantities 
can be measured from afterglow data. We thus have an opportunity 
to check if the values of cr obtained from observations of GRB af- 
terglows are consistent with the y-ray efficiencies £y measured for 
the same GRBs. Carrying out this test is our goal. 

In ^we write down the standard jump conditions for a mag- 
netized relativistic "perpendicular" shock in which the magnetic 
field is perpendicular to the flow velocity (or parallel to the shock 
front). By solving the jump conditions, we calculate the efficiency 
with which the kinetic energy of a cold relativistic magnetized fluid 
is converted by the shock to thermal energy. In ^ we calculate the 
efficiency of the internal shock model and compare it against ob- 
servations, and in ^ we carry out a similar exercise for the reverse 
shock model. In both cases, we show that there is an inconsistency 
between the predictions of the model and measured values of £y 
and cr. We discuss the implications of this result in ^and suggest 
possible solutions. 



2 RELATIVISTIC PERPENDICULAR SHOCK 
2.1 Preliminaries 

The problem o f interest was discuss ed in detail by 
iKennel & Coronitil jl984 hereafter lKC84h . We follow their 
methods with a few minor changes. We consider a cold magnetized 
fluid with a magnetic field strength Bo in its rest frame. We assume 
ideal magnetohydrodynamics (MHD) and set the electric field 
in the rest frame to zero. Transforming to a frame in which the 
magnetized fluid moves with dimensionless velocity p = v/c, 
Lorentz factor y = 1 / -^1 in a direction perpendicular to the 
magnetic field, the magnetic and electric fields become 



B = jBo 



uBq = -B, 
7 



:/?y = (y2- 1) 



1/2 



(1) 



where u is the relativistic 4-velocity. The fields B and E in the new 
frame are parallel and perpendicular, respectively, to So, and each 
is also perpendicular to the velocity. 



' It would be very difficult to produce the observed y-ray variability in the 
rever se shock model unless there is relativistic turbulence in the shocked 
fluid iNaravan & Kumail200g| ; lLazar. Nakar & Piranl200^ . 



We define the magnetization parameter cr of the moving fluid 
as the ratio of the Poynting energy flux to the particle rest energy 
flux. Thus 



cEB/4n B^IAn Bl/4n 



(2) 



where n is the particle number density in the fluid rest frame and 
m is the mass of each particle. In the final expression, the numer- 
ator is the rest frame "enthalpy" of the magnetic field, which is 
equal to ITb/(Tb - l)](BQ/Sn) (taking the adiabatic index of the 
magnetic field = 2 for compression transverse to the field), 
and the denominator is the rest energy density. Since we see that 
cr depends only on rest frame quantities, it is a relativistic invariant 
and is frame-independent. KC84 use a slightly different definition 
of cr where they replace y- in the denominator of the third quan- 
tity in equation ^ by yu. As a result, their cr is not truly frame- 
independent. However, the difference between the two definitions 
is negligibly small for highly relativistic flows. 



2.2 Jump Conditions 

We follow !kc84 except that we use the definition of cr given in 
equation l|2) and avoid certain approximations. We use subscript u 
for the gas upstream of the shock and subscript d for the down- 
stream gas. The upstream gas is cold (P,, = 0), has a magnetization 
parameter cr, rest number density n„, and moves with Lorentz factor 
y„ in the frame of the shock. The downstream gas is hot (Pj t 0) 
with adiabatic index F, has number density nj, and Lorentz factor 
yj. In the shock frame, the magnetic fields in the two regions, 6„ 
and Bj, are related by 



B = He = —E = — s 
I'd lid y„iij 



(3) 



where we have used the fact that the electric field is continuous 
across the shock. 

The upstream and downstream gas enthalpy per particle are, 
respectively. 



1 + 



(F — 1 ) rijmc^ 



(4) 



The second term inside the square brackets is a dimensionless num- 
ber which describes the thermal enthalpy per particle of the shocked 
gas. It can be written as 



(F - 1) njmc^ 



kTd_ 



T- 1 

(5) 



where O^i is the relativistic temperature of the downstream gas. 
When Oil <k 1, the gas is non-relativistic, and we have F = 5/3, 
h{6d) = 5/2, whereas when dj » 1, the gas is ultra-relativistic, 
and we have F = 4/3, h(6d) = 4. At intermediate temperatures 
(Oil ~ 1), h(6j) can be written in terms of modified Bessel functions 

I w 1 

(see Chandrasekhar 1960). For simplicity, we use the following ap- 
proximation (Service. 1986,) , 

10 -H 200 



h(9). 



4 + 56 



(6) 



which is sufficiently accurate for our purposes. In principle, if the 
jet material consists of a normal electron-proton plasma, we should 
allow for two species of particles in the shocked gas, each with a 
different temperature. We ignore this complication for simplicity. 
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We have three jump conditions across the shock, correspond- 
ing to three fundamental conservation laws. First, mass conserva- 
tion implies that the mass fluxes on the two sides of the shock must 
be equal, i.e., 

n„M„ = njiirf. (7) 
Energy conservation requires the energy fluxes to be equal, i.e.. 



EB„ 



■ ndUdJdt^d + 



EBd 



An An 
Finally, momentum conservation gives the condition 



2 ^(Bt + Et) 



: ndUj^id + Pd + 



(8) 



(9) 



In the last equation, the terms involving the electric field cancel 
since £"„ = Ed- Eliminating jid between equations {8]( and ^ and 
simplifying, we obtain the following expression for 0^: 



ft 



7i/ 



(10) 



y'd"" 

Ud Jd! \ 2u„ 2ud ) 

In addition, equation ([8} can be rewritten in the following simplified 
form. 



1 + h(ed)ed = —(1 + cr) - — cr. 

yd lid 



(11) 



Given the upstream quantities u„ and cr, it is straightforward to 
solve equations jlOj and i ll II ) numerically. We guess a value for the 
downstream velocity Ud and calculate 6d using equation jlOj . We 
then compute h(dd) using the approximation (|6j and check whether 
the condition jilt is satisfied. If it is not, we numerically adjust Ud 
until the condition is satisfied. We then have the complete solution 
for all downstream quantities: jd, Ud, "d/n,,, 6d, /id/mc^, Bd/B,,. 

The results presented in the following sections use the above 
numerical approach to solve the jump conditions. An alternate ap- 
proach is to make suitable approximations and obtain analytical 
solutions of the jump conditions. Appendix lAl presents analytical 
solutions corresponding to a number of useful limits. 



3 INTERNAL SHOCK MODEL 
3.1 Solving the Jump Conditions 

We consider two identical blobs, each with magnetization cr, ap- 
proaching each other and colliding. In the center of mass frame, 
the blobs have Lorentz factors y and relativistic velocities u = 
- Vt" ~ 1- As a result of the collision, two identical shocks move 
(in opposite directions) into the two blobs. 

For given values of y and cr, we solve the jump conditions nu- 
merically and calculate all quantities of interest in the shocked gas. 
We begin by assuming a value for the upstream Lorentz factor y„ in 
the frame of one of the shocks. Following the procedure described 
in 32.21 we solve for the downstream Lorentz factor yd- From y,, 
and yd, we calculate the relative Lorentz factor y,,^ between the 
two regions (relativistic velocity subtraction) and check whether it 
corresponds to the desired value of y. If not, we adjust y„ until we 
obtain y„<; - y. We then have the solution. 

Having obtained the solution, we switch to the rest frame of 
the shocked gas. We assume that a fraction 6^ of the thermal en- 
thalpy of the shocked gas Wg.^s goes into electrons and that it is 
entirely radiated in y-ray^ This gives the energy Ey that goes into 



This is perhaps a little optimistic. It is possible that only the gas internal 
energy is radiated, which is 1 /F times the enthalpy. 



-1.5 




-1.5 




Figure 1. Left: y-ray efficiency 6y vs the upstream magnetization pai'ame- 
ter o" for an internal shock between two identical biobs. It is assumed that 
a fraction = 0.2 of the enthalpy in the shocked gas is radiated as prompt 
y-rays. From below, the curves correspond to blob Lorentz factor in the 
center-of-mass frame of y = 1.05, 1.1, 1.2, 1.4, 2, 4, 10, respectively. Equiv- 
alently, the inter-blob Lorentz factor between the two blobs is yib = 1.21, 
1.42, 1.88,2.92,7,31, 199, respectively. The symbols refer to obsei'v ational 
data. The filled circles con'espond to the first nine GRBs in Table[T]/PK0^ 
sample) and the open circles to the last three GRBs. The arrows on the sym- 
bols indicate that the estimated values of cr are lower limits (eq.[T8). Right: 
CoiTesponding results for the case when = 1, i.e., all the thermal enthalpy 
in the shocked gas comes out in prompt y-rays. 



y-rays. The remaining unradiated enthalpy, which consists of rest 
mass enthalpy VKcst, remaining gas thermal enthalpy (I - £e)W^as 
and magnetic enthalpy Wg, contributes to the kinetic energy Eq that 
goes into the afterglow. Thus, we estimate the efficiency £y of y-ray 
emission, the fraction of the total energy that goes into y-rays, to 
be 

Ey QW„„, 
£y = = . (12) 

" Ey+Eo Wgas + M/,es, + IVs 

The quantities W^^s and Wy^^t are easily obtained from the 
shock solution. Per particle, they are given by 



■ h{0d)9dmc~ , Wrest = 



(13) 



To calculate Wg, we first need to calculate the magnetic field of 
the shocked gas in the center of mass (CM) frame of the colliding 
blobs. This is given by 



Bd 

BcM - — 



-s„. 



yd yi,Ud 
Then the magnetic enthalpy per particle is 



Wb 



BIm 
Antid 



I'm 2 

— crmc . 

Ud 



Thus we obtain 

€MOd)ed 



(14) 



(15) 



(16) 



l+h(ed)Od + (Uu/ud)cr' 

Although all quantities have been estimated in the rest frame of 
the shocked gas, it is easily shown that Lorentz transforming to a 
different frame, e.g., the observer frame, will leave Sy unchanged. 

For a given value of e^, the y-ray efficiency depends on two 
parameters, the magnetization cr of the blobs and their Lorentz fac- 
tor y in the center-of-mass frame. Figure [T] shows how varies 
as a function of cr for selected values of y, as listed in the figure 
caption. Instead of y, we could express the results in terms of the 
relative inter-blob Lorentz factor of the blobs yib. The values of yib 
are also given in the figure caption. 

Figure[T]indicates that the maximum radiative efficiency is ob- 
tained for unmagnetized blobs. As the magnetization increases, the 
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amount of thermal energy generated in the shock decreases, reduc- 
ing th e radiative efficiency. This result was already discussed by 
IkC84| . The analytical approximations in Appendix |A] give further 
details. From 9 A3 1 we see that, for cr <c 1, the enthalpy of the 
shocked gas (which is proportional to dj) varies as 1 - (5 /2)ct, i.e., it 
reduces with increasing magnetization. The reduction is quite pro- 
nounced once (T > 1; for cr » 1, the enthalpy scales as l/cr QAlt . 

Note, however, that there is always a shock solution for any 
choice of y and cr. This may appear a little surprising since, as we 
show in Appendix [A] a magnetized shock is possible only if the 
upstream velocity m„ exceeds ^/lT. Thus, for strongly magnetized 
blobs moving in the center-of-mass frame with relatively low ve- 
locities, one might think there should be no shock. For example, 
for 7 = 10 or u = 9.95, one might expect the shock to disappear 
once cr > 3.15, whereas Fig. [T] shows results for cr values much 
above this limit. 

The explanation is simple. What is relevant for the existence or 
otherwise of a shock is the upstream velocity in the shock frame, not 
the center-of-mass frame. In all our solutions, the shock moves out- 
ward in the center-of-mass frame. When » cr, the shock veloc- 
ity is not very large. However, once cr exceeds ir, the shock moves 
quite rapidly into the blob. In fact, it moves so rapidly that the ve- 
locity of the upstream gas k„ as seen in the shock frame becomes 
larger than ^fa, thus permitting a shock. Shocks in this regime are 
however weak (see 3A2b . and their radiative efficiencies Ey are cor- 
respondingly very low. As an example, note in Fig.[T]the very low 
efficiency of the y = 10 model (the uppermost curves in the two 
panels) when cr = 10. 

3.2 Comparison with GRB Data 

The data we use are listed in Table [Ti fPanaitescu & Kum^ j2002l 
hereafter PK02) have analyzed afterglow observations of ten GRBs, 
and have derived for these objects the parameters we need. We 
make use of the results in Tables 2 and 3 of their paper. We include 
GRB 970508 for which we take from lBloom. Frail & Kulkamil 
( l2003h , but we omit GRB 980159 since the red shift is not known. 
For the remaining 8 systems, we take |PK02l 's estimates of Ey, 
the beaming-corrected y-ray emission in the 20-2000keV band, 
and Eq, the beaming-corrected kinetic energy of the afterglow, 
and compute the y-ray efficiency parameter 6y. The afterglow data 
considered by PK02 did not include observations during the early 
stages of the afterglow (first day or so). During this time, the exter- 
nal shock is expected to be somewhat radiative. It is thus possible 
that PK02 slightly underestimated Eq in their models. To allow for 
this, we double their values of Eq and estimate the y-ray efficiency 
by Cy ~ EyIiEy + 2Eo). (The correction factor of 2 is probably too 
large, but our aim is to be conservative.) 

In addition, we also estimate the magnetization cr of the 
jet ejecta. Based on numerical a nd analytical work on cold 
magnetically-accelerated GRB jets, iTchekhovskov et al] ( l2010bl) 
have shown that the following inequality must be satisfied, 

y^ sinSj < 15o-"^ (17) 

where yj and 9j are the Lorentz factor and opening angle of the 
jet ejecta at the conclusion of the prompt emission phase, i.e., just 
before the onset of the a fterglow phase. The factor of 15 is a loga- 
rithmic term. Since |PK02I have estimated yj and Oj for the 9 GRBs 
of interest to us (see Table [T}, from their data we obtain for each 
GRB a lower limit on cr. 



The solid circles in Fig. \T \ show the values of cr„,|„ and £y for the 
nine GRBs from the |PK02l sample. 

In addition, we have gone through the literature and esti- 

GRB 021004 (relevant 
Li & Chevalier] 
2009 



mated £y and cr,^in for three more bursts 



data taken from iBloom. F rail & Kulkarni 2003; 
2003), GRB 080916C IXbdo et al. 2009; Greine r et al. 
Kumar & Bamiol Duran 2009), GRB 090510 l lAckermann et al.1 
l2010i ; iKumar & Bamiol Duranii201ol) . These three GRBs are shown 
in Fig. [T] with open circles, and the corresponding data are given in 
the last three lines of Table[T] 

In order to compare the data with the predictions of our in- 
ternal shock model, we need to assume a value for the electron 
heating parameter e^. From modeling afterglow observations it is 
possi ble to estimate for the f orward shock in individual GRBs 
(e.g.. Pana itescu & Kumaill2002h . The median value from a sample 
of 39 GRB afterglows is ~ 0.2 (Santana & Bamiol Duran, in 
preparation). The theoretical curves in the left panel of Fig.[T]cor- 
respond to this value of e^. A quick look shows that the model com- 
pletely misses the observations for all twelve GRBs in our sample. 
It has been known for some time that internal shocks involving un- 
magnetized shocks can not easily achieve the y-ray efficiencies re- 
quired by observations ( Kumailll999l; lLazzati. Ghisellini & Celottil 

-oili 



1 7j sin Sj 
15 



(18) 



I I999I ; ^Panaitesc u. Spada & Meszarosll l999l) . Magnetization makes 
the problem worse. For the cr values estimated for GRBs, the pre- 
dicted efficiency is much lower than for the unmagnetized case, 
so the discrepancy is quite large. Note further that the values of cr 
plotted for the individual GRBs correspond to cr„,i„ (eq.[T8}. This 
means that the points might actually lie even farther to the right, 
which would make the discrepancy impossibly large. 

The panel on the right in Fig. [T] shows the highly optimistic 
case when q = 1. This might coiTespond, for example, to an 
electron-positron jet. Even in this case, the majority of GRBs are in- 
consistent with the magnetized internal shock model. We thus con- 
clude that, if jet ejections are magnetized and are described by ideal 
MHD, and if the blob Lorentz factors y and yit are not very differ- 
ent from the values considered in Fig.[T|(note that the largest value 
considered is yib ~ 200 which is very unlikely to be exceeded), 
then the magnetized internal shock model is ruled out conclusively 
by the data. 



4 REVERSE SHOCK MODEL 
4.1 Solving the Jump Conditions 

In the case of t he reverse shock, w e have to consider four regions, 
as described in lSari & PirMl ( ll995L hereafter |SP95|) : 

1. The external unshocked ISM, which is at rest in the lab frame 

2. The shocked ISM 

3. The shocked jet ejecta 

4. The unshocked jet ejecta 

Following ISP95L we use subscripts I, 2, 3, 4 to identify quantities 
in the four regions. Regions 2 and 3 are in pressure equilibrium 
across the contact discontinuity and move with the same Lorentz 
factor. As measured in the lab frame, region 1 is at rest, regions 
2 and 3 move with Lorentz factor y^ (c for contact), and region 4 
moves with the jet Lorentz factor yj. 

In the analysis below we consider several distinct frames. 
First, we have the lab frame in which the unshocked jet moves with 
Lorentz factor y^ and the shocked gas moves with y^. Next, we have 
the frame of the forward shock in which the unshocked and shocked 
ISM move with Lorentz factors yi and y2. Then, we have the frame 
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Table 1. GRB data (see ^S^Jfor details). 



Source 








Ji 






Ref. 


GRB 970508 


3.8 


20. 


0.087 


150 


18.3 


9.86 


1,2 


GRB 990123 


4.9 


1.5 


0.62 


300 


2.1 


0.54 


1 


GRB 990510 


1.3 


1.4 


0.32 


140 


3.1 


0.25 


1 


GRB 991208 


18 


2.4 


0.79 


68 


12.8 


1.01 


1 


GRB 991216 


3.0 


1.1 


0.58 


150 


2.7 


0.22 


1 


GRB 000301c 


6.6 


3.3 


0.50 


160 


13.7 


6.38 


1 


GRB 000418 


148 


32 


0.70 


90 


50.0 


21.1 


1 


GRB 000926 


15 


3.2 


0.70 


130 


8.1 


1.49 


1 


GRB 010222 


11 


5.1 


0.52 


110 


4.6 


0.35 


1 


GRB 021004 


560 


400 


0.58 


55 


12.7 


0.65 


2,3 


GRB 0809 16C 


8.8 X 10'' 


105 


0.47 


> 880 


> 2.2 


> 5.07 


4,5,6 


GRB 090510 


1100 


4000 


0.22 


> 1200 


0.7 


> 0.96 


7,8 



" Beaming-coiTected energy in y-rays in units of 10 erg 

* Beaming-corrected kinetic energy in the afterglow in units of 10^" erg 

' Calculated using eq. \\2\ for the last three GRBs, but with Eq replaced by 2£o for the 

first nine GRBs 

R eferences: 1 |PK02l: 2 ^iBloom. Frail & Kulkaiir] l2003t): 3 -lU & Chevalieil j2003l) : 4 
-Ubdo et al. 1 2009|); 5 ^ [Grein er et al. ( 2009): 6 Kumar & Bamiol Duraj j20O9[) ; 7 - 
lAckermann et al., t2010t) : 8 - iKumar & Barniol Duran c2Q10i) . 



of the reverse shock in which the unshocked and shocked jet ejecta 
move with Lorentz factors 74 and 73 . Finally, we have the frame of 
the shocked gas. This frame moves with respect to the lab frame 
with a Lorentz factor y^., with respect to the forward shock frame 
with a relative Lorentz factor 72, and with respect to the reverse 
shock frame with a relative Lorentz factor 73. 

We consider regions 1 and 2 to be essentially unmagnetized 
and treat the forward shock between these two regions as a hydro- 
dynamic shock. Furthermore, we assume that the relative Lorentz 
factor across this shock is large, i.e., y^, m^. » L Let us trans- 
form into the frame of the forward shock. The upstream gas is cold 
{P\ = 0, recall that subscript 1 refers to region 1), is unmagnetized 
(by assumption, cr = 0) and has a large Lorentz factor (71 » 1). 
We can use the results given in gA3l (with subscripts u, d replaced 
by 1 and 2) in the limit cr —> to calculate the properties of the 
shocked gas. We then obtain the following standard results for the 
downstream gas, 



1 

3' 



_3 

2V2 



— =, n2 = 2V2 



P2 = -y]rnmc\ (19) 



The relative Lorentz factor between the two regions, which we call 
jc, is equal to 71 / V2. Thus, we find 

4 



"2 = 47t.ni, 



(20) 



Consider now the reverse shock between the magnetized re- 
gions 4 and 3. The jet ejecta have a Lorentz factor yj in the lab 
frame and a magnetization cr. In the frame of the reverse shock, we 
do not know a priory the value of the upstream Lorentz factor 74. 
Therefore, as in ^ we will solve for 74 via the jump conditions 
(all the relations given in ^are valid, except that subscripts u and 
d should be replaced by 4 and 3, respectively), plus an additional 
requirement. In ^ the additional constraint was the value of 7 (or 
equivalently 711,). Here, for easy comparison with previous work in 
the literature, we will fit a target value of the "relativity" parameter 
^ defined injSP95 and Giannios et al. (2008). This parameter is less 
than unity for a relativistic shock and greater than unity for a New- 
tonian shock. In Appendix |B] we show that for a magnetized flow 



3 114 
y] m 



(l+o-) 



1/2 



(21) 



where Ri^^ and R, are the deceleration radius and the spreading 
radius of the expanding ejecta. 

Given the jet Lorentz factor yj in the lab frame, the magnetiza- 
tion of the jet material cr, and a target value of the relativity param- 
eter ^, the calculation proceeds as follows. We begin by guessing a 
value for 74, the Lorentz factor of the upstream jet ejecta as viewed 
in the frame of the reverse shock. Then, as described in ^ we solve 
for all quantities in the downstream region 3. In the rest frame of 
the shocked gas, the gas pressure is equal to 



^3 gas = njO^mc^ = n^ — O^mc^, 

and the magnetic pressure is equal to 



3,mag 



8;r 



"4 -^4 "4 7 

— = ^4 — -crmc 



(22) 



(23) 



Since regions 2 and 3 are in pressure balance, we thus obtain the 
following condition. 



: «4 



luj 



(24) 



The Lorentz factor on the right-hand side of equation ( I24t is 
straightforward to calculate. In the reverse shock frame, we know 
that region 4 has a Lorentz factor 74, whereas its Lorentz factor 
in the lab frame is yj (which is given). Thus, we can calculate the 
Lorentz factor of the reverse shock as seen in the lab frame (rel- 
ativistic velocity subtraction). Once we have this quantity, we can 
transform 73 (which is in the reverse shock frame) to the lab frame 
to obtain y^, the Lorentz factor of region 3 (as well as region 2) in 
the lab frame. 

Having calculated y,., we obtain the density ratio n^jn^ from 
equation I l24t . and hence the value of ^ from equation ( I21t . We 
check this against the target value of ^, and numerically adjust 74 
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Figure 2. Similar to Fig.[T] but for a magnetized reverse sliock. The curves 
correspond to jet Lorentz factor -yj = 300 as measured in the lab frame, and 
(from below) relativity parameter f = 3.162, 1, 0.3162, 0.1, 0.03162, 0.01. 
The panel on the left is for £g =0.2 and that on the right for = 1 . 



until we achieve the value of f we seek. At this point, we have the 
solution to the problem. 

Once we have the solution, we can calculate the parameter 6y. 
As before, we will assume that a fraction 6^ of the gas thermal en- 
thalpy in region 3 is radiated. As measured in the lab frame, this 
corresponds to an energy per particle of e^yc /i(93)03mc-^. To cal- 
culate the total energy per particle of the system, it is simplest to 
consider the unshocked jet fluid. In its own rest frame, the enthalpy 
per particle is mc^(l + cr), and this gets multiplied by yj when we 
transform to the lab framed Thus we obtain 



yj(l+iT) ■ 



(25) 



Apart from e^, this reverse shock model has three parameters: 
the jet Lorentz factor yj, the magnetization parameter of the jet ma- 
terial cr, and the relativity parameter ^. Figure |2] shows results for 
a fixed value of = 300 (the results hardly change for other val- 
ues, e.g., 100 or 1000). The curves correspond to selected values 
of the relativi ty parameter As in the c ase of unmagnetized re- 
verse shocks ( ISP95l iGiannios et al.ll2008l) , we find that the y-ray 
efficiency is highest for highly relativistic shocks ^ 1), and the 
efficiency is very poor for Newtonian shocks (f > 1). In addition, 
for a given value of ^, the efficiency decreases as the magnetiza- 
ti on increases. This result is in qualitative agree men t with the work 
of Zhang & Kobavashj j2005l) . lLvutik"ovl j2005l) and lGiannios et all 
( l2008l) . 



4.2 Comparison with GRB Data 

As in the case of the internal shock model, a comparison of the 
predictions of the magnetized reverse shock model with GRB data 
(Fig.[2j indicates that the model has no hope of satisfying the obser- 
vations. For £f = 0.2, which we consider a reasonable value, not a 
single GRB agrees with the model even if we assume an extremely 
relativistic shock with ^ = 0.01. For = I, which in our opinion is 
rather optimistic, a few systems do fall inside the model curves, but 
far too many systems still remain unexplained. We thus conclude 
that, if jet ejections are cold and magnetized and are described by 
ideal MHD, then the magnetized reverse shock model considered 
here is ruled out by the data. 



5 SUMMARY AND DISCUSSION 

The magnetic acceleration paradigm for relativisti c jets i s theoret- 



i97e 


; 


2001 






reference list). Recent advances in numerical techniques, cou- 
pled with analytical methods, have led to a deeper understand- 
of how Poynting-dominated jets accelerate to l arge Lorentz 



ing ot now Jr'oynting-aommatea lets accelerate to l arge Lorentz 
factors jKoinissarov l2004l: iKomi ss aroy et all [20071 12009|, |201(]|: 



iTchekhovskov et al 



2008l ]2009, 2010aHb|). In the specific context 



of ultra-relativistic GRB jets, it has recently been shown that the 
collimation angle Oj and the Lorentz factor yj of a magnetically 
accelerated jet are no t independent but are relate d via the magne- 
tization parameter cr jTchekhovskov et al]|20 1 Obi) . Estimates of 6j 
and yj of GRB jets, obtained by modeling afterglow data, are gen- 
erally consistent with these jets having cr ~ 1 just prior to the onset 
of the afterglow jTchekhovskov et al]|2010b.) . This indicates that 
GRB jets successfully convert about half of their initial Poynting 
flux to matter kinetic energy by the time they reach the decelera- 
tion radius. These jets are thus energetically efficient. 

Whe n it com es to radiative efficiency, however, cr ~ 1 is not 
sufficient. IKC84 showed that perpendicular shocks in cold mag- 
netized gas produce thermal energy very inefficiently unless cr is 
much less than unity. We have explored this issue in detail in the 
context of GRB internal shocks and reverse shocks. Both of these 
shocks occur in the material ejected in a GRB and are expected to 
be magnetized. The geometry is also such that the shocks will be 
perpendicular, i.e., the magnetic field will be perpendicular to the 
velocity vector, or parallel to the shock front. We have analyzed 
such shocks, assuming that a fraction of the gas thermal energy 
in the shocked gas goes into electrons and that this energy is en- 
tirely radiated in prompt radiation. We consider two values of 6^, 
viz., 6^ = 0.2, which we consider to be a reasonable estimate, and 
6p = 1, which is highly optimistic. 

Our calculations indicate that, once cr exceeds about 0.1, the 
efficiency of thermalization begins to fall noticeably, and that the 
drop becomes quite precipitous once cr > 1 (Figs.[T][2j- GRB ob- 
servations indicate that the prompt y-ray emission is quite efficient, 
with the efficiency parameter 6^ (defined in eq. [T2t being typically 
of order 0.5 or larger (Table[T]. On the other hand, not a single GRB 
has cr < 0.1, as needed to obtain such high efficiency in a cold mag- 
netized shock, and half our sample has cr > 1, where radiative effi- 
ciency is very poor. The implication is that GRB prompt emission 
cannot be produced by either internal shocks or the reverse shock, 
if jets are cold and magnetically accelerated. This conclusion is 
hard to avoid. Even with very optimistic assumptions, e.g., all the 
thermal energy goes into electrons (e^ = 1, which might happen if 
the jet is made entirely of electrons and positrons), and is immedi- 
ately radiated in y-rays, the calculated efficiency is far below what 
is needed to explain the observations. 

We consider here several possible resolutions of this puzzle, 
none of which is very satisfactory^ 

One possibility is to associate the prompt emission with the 
forward shock, which is very weakly magnetized (cr « 1) and 
therefore converts a large fraction of the jet kinetic energy into ther- 
mal energy. This alone is not enough since the thermal energy must 



^ Note that, for both of the above quantities, we first calculate the enthalpy 
in the rest frame of the gas, where there is no net momentum. Thus, trans- 
formation of the energy to another frame requires only multiplication by the 
relevant Lorentz factor 



'* IZhang & Yai] j201lh have considered internal collisions for a magneti- 
cally dominated outflow and suggest that this could facilitate dissipation of 
magnetic energy via reconnection. We do not discuss this particular process 
here. 
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Figure 3. Similar to Figs.[T]and|2] except it is assumed tliat, in addition to a 
fraction = 0.2 of the gas thermal enthalpy, the entire magnetic enthalpy is 
also radiated as prompt y-rays. The panel on the left is for an intemal shock 
and that on the right for a reverse shock. 



then be radiated with greater than 50% efficiency in order to ex- 
plain the observed values of €y. Therefore, nearly all the thermal 
energy should go into electrons (e^ ~ 1), which is very unlikely 
for the electron-proton plasma we expect to be present in the for- 
ward shock. In any case, the forwar d shock has been con vincingly 
associated with afterglow emis sion dOehrels eT al. 2009; MeszarosI 
I2OO2I : |Piranl2005l : IZhang|l2007l : and references therein), and it does 
not seem likely that the same region will also produce the prompt 
emission. 

A second possibility is that much of the magnetic energy in 
a GRB shock is somehow converted to particle thermal energy. 
That is, when cr is large and most of the energy density in the post- 
shock gas is in the form of magnetic energy, there is a mechanism 
whereby this energy is converted to particle energy. A scenario 
where this can happen is if the pre-shock gas is "striped" as in cur- 
rent models of ener gy dissipation in the magn etized wind of pulsars 
jLvubarskv & KirkllS)Oll ; lLvubarskvll2010al). A striped morphol- 
ogy i s not obvious for a GRB jet (but see McKinnev & Uzdenskvl 
I2OIOI) . However, if it is present, we do expect a substantial frac- 
tion of the magnetic energy to be dissipated in the shock. Figure |3] 
shows results for a hypothetical model in which we assume that, in 
addition to a fraction = 0.2 of the gas thermal enthalpy, 100% of 
the magnetic energy is radiated. This model does explain the GRB 
observations but at the price of making a very extreme (and theo- 
retically unsupported) assumption. We do not endorse this model 
but present it merely as a way to emphasize how difficult it is to 
explain the radiative efficiency of GRB prompt emission. 

A third possibility is that the magnetic energy is dissipated, not 
through a shock, but through some other "current-driven" mecha- 
nism such as reconnection. Poynting-dominated magnetically ac- 
celerated jets are fairly stable once they are ultra-relativistic (e.g., 
iNaravan. Li & Tchekhovskovll2009l) and are unlikely to have large- 
amplitude fluctuations that might drive reconnection. However, it 
is conceivabl e that these jets lose their stab ility once they reach a 
large radius dMcKirmev & Uzdenskvl I20 id) , e.g., the deceleration 
radius where the jet meets the external medium and begins to slow 
down. Whether the instability would be powerful enough to drive 
wholesale reconnection and convert most of the magnetic energy 
into particle energy is an open question. As Fig. |3] shows, some- 
thing like this is needed if one is to explain the data. 

Another possibility is that our assumption of cold gas, 
whose acceleration is entirely by magnetic means, is in- 
correct. N on-relativistic MHD simulations of m agnetized jets 
( iMoll. Spr uit & Obergaulinger 2008; Molil |2009|) indicate that 
these jets develop a kink instability which might lead to dissipa- 
tion. We could then have a scenario in which the jet starts off mag- 



netically dominated at the base but quickly dissipates its magnetic 
energy into heat while the jet is still non- or quasi-relativistic. Fur- 
ther acceleration of the jet is then driven by the thermal pressure 
of the heated gas. Thus, we no longer have a magnetically driven 
jet, but something akin to the standard fireball model of a GRB. 
Clearly, the calculations presented here, which are restricted to cold 
magnetized gas, are not relevant for such a model. 

Finally, it is possible that the prompt emission in GRBs 
is not produced in the jet at a large distance from the progen- 
itor, but rather in the photospheric region where the jet ejecta 
first become transparent. Models of this form have been devel- 
oped fThompsorJI 19941; iMeszaros & Reesll2000l; I Rees & MeszarosI 



I2OO5'; 'Pe'er. Meszaros. Reesl l2006l; iGianniosI I2OO8I ; iBeloborodov 
2010; Metzger et al. 201 1() and it is claimed that they prO' 
duce prompt y-ray emission with high radiative efficienc; 
and with the correct spectr um 



IVurm. Beloborodov & Poutanei 



1 mgn radiative en i ciency 

(e.g., IPe'er & Rvdel bOlll 

20 111). Magnetic fields may play 



a role in photospheric models t Uzdenskv & McKinnev luOlO h. but 
the role of shocks is unclear. Our analysis is not applicable to these 
models. 
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APPENDIX A: ANALYTICAL APPROXIMATIONS FOR A 
PERPENDICULAR SHOCK 

A cold hydrodynamic flow can have a shock for any choice of the 
upstream velocity u„. This is because no signals can propagate in 
the cold gas and so the upstream gas is always supersonic. A mag- 
netized fluid is diff'erent. Even if the gas is cold, Alfven and fast 
magnetosonic waves can still propagate in the fluid. Thus a shock 
is possible only if the upstream gas moves faster than these waves. 

For the particular geometry we have considered, viz., a per- 
pendicular shock with magnetic field perpendicular to the veloc- 
ity vector, the relevant wave speed is that of the fast magnetosonic 
wave, which is given (in the comoving frame of the gas) by 



^fms 



o- + 1 "™ 

Thus we can have a shock only if 



(Al) 



(A2) 



We now consider a number of limiting cases. When iil » 
o", we expect to have a strong shock, whereas when ul is only 
marginally greater than cr, we expect a weak shock. In addition, we 
have different results depending on whether » I (ultrarelativis- 
tic) or iq <c 1 (nonrelativistic), and on whet her cr » 1 (strongly 
magnetized) or tr « 1 (weakly magnetized). IKC84| considered a 
couple of important cases, but here we present scalings for all the 
different regimes. Figure lATI identifies the regimes and labels them 
by the respective subsection where each is discussed. 
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Figure Al. Shows different regimes for a perpendicular shock in cold up- 
stream gas. The horizontal axis indicates the magnetization parameter o" of 
the upstream gas (eq.|2] and the vertical axis indicates the square of the up- 
stream relativistic velocity m„ as measured in the frame of the shock. There 
is no shock when w„ < -\/(T since the upstream gas moves slower than 
the fast magnetosonic wave speed Uf^j = Vcr (eqs. Al, A2). The regions 
marked A.1.,...,A.5, in the plot correspond to different shock regimes. The 
labels refer to the section numbers in Appendix IaI where each is discussed. 
Regions A.l. and A. 2. correspond to highly magnetized shocks (cr 2> 1), 
which are necessarily also highly relativistic (u„ 2> 1). A.l. describes a 
strong shock, where u„ 2> Vcr, while A. 2. describes a weak shock where 
Mu is only marginally greater than V^?- Regions A. 3., A. 4., A. 5., correspond 
to weakly magnetized shocks (cr « 1). A. 3. is highly relativistic (u„ 2> 1), 
where the shock is necessarily strong. A. 4. is non-relativistic (u„ <K 1), but 
still with a strong shock, while A. 5. corresponds to a non-relativistic weak 
shock, where h„ is only marginally greater than 



A2 Weak Relativistic Shock with Strong Magnetization: 

K, ^ o- » 1 

As approaches cr, the shock becomes progressively weaker. Let 
us write = cr(l 4- A), with A <c 1. As A becomes progressively 
smaller, less and less of the kinetic energy of the upstream gas is 
thermalized in the shock. In this limit, the downstream temperature 
6,/ becomes non-relativistic and so we set hidS) = 5/2. In this limit, 
the leading terms in the solution are as follows: 



A 44 , 
1-- + — 1 

3 81 



fi„ 30-' 



3 

Slcr" 



-A' + o 



4 1 1 \ , 8 / 4 \ , /A^ An 
~ — 1 - - + -7 1 + o\—-—\ 

81 I o- cr'- 8ll 3cr Icr^' cr ) 



(A6) 
(A7) 
(A8) 



Note that the leading term in the temperature of the shocked gas 
goes as a', i.e., the shock is extremely inefficient. This is a charac- 
teristic feature of weak shocks. 



A3 Strong Relativistic Shock with Weak Magnetization: 

ir » 1 » cr 



This case has again been considered bv lKC84 With h(ea) = 4, the 
results are 



,19 9 2 
'ij X — I — cr cr + 



o{cr^). 



— « 3 - llcr + 96cr^ + o{cr-) , 

Sd 1 /, 5 111 ,\ / 2\ 

— a — 1 cr H cr- M- o cr^ . 

"„ 3V2\ 2 S I ^ > 



(A9) 
(AlO) 
(All) 



Again, our result for dj differs from that given in lKC84l The results 
for an unmagnetized relativistic shock are recovered by simply set- 
ting cr = in the above relations. 



Al Strong Relativistic Shock with Strong Magnetization: 

m: » cr » 1 



This case has been considered by IkC84| who give the following 
results: 



1 



1 



M , a; cr 4 1 h o — 

8 64cr Icr 

5rf , 1 3 

— a 1 + + o\ 

B„ 2cr 16cr2 



u„ 



1 



3 
16cr 



r3/2 



(A3) 
(A4) 
(A5) 



where o{x) denotes terms of higher order than x, i.e., o{x) /x ^ 
as X — » 0. In the above, we remind the reader that the magnetic 
field strengths S„ and are measured in the shock frame. These 
results are obtained by expanding the jump conditions as power se- 
ries in the small quantity 1/cr, and matching terms of similar order. 
We have set hiOj) = 4, as appropriate for relativistically hot dow n- 
stream gas. Our result for Oj/u,, differs from that given in lKC84 



A4 Strong Non-Relativistic Shock with Weak Magnetization: 

l»l3l»cr 

Here we consider the non-relativistic case and replace m„, uj by 
Pd- Also, we set HOd) = 5/2. Then we find 



9 cr 9cr2 



^.4-18^ + 117^ + 44 

Bu PI Pt \Pt 



Od 



16 



■A 



21 
16 



cr + o(cr) . 



(A12) 
(A 13) 
(A14) 



The solution for an unmagnetized shock is obtained by setting cr = 
0. As an aside, we note that cr is related to the Alfven wave speed 
va by 



Ammc^ 



-Pi 



(A 15) 
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AS Weak Non-Relativistic Shock with Weak Magnetization: 

Finally, we consider the case when the shock is non-relativistic and 
j8,, = Vo^(l + A) with y/(T <c A <c 1. In this limit, we find 



— » 1 + -A A^ + o A^ , 

B„ 3 81 ^ ) 

0„.pA3_ii£A^]^ + o(AV). 
\81 243 / ^ ' 



APPENDIX B: THE RELATIVITY PARAMETER f 



(A16) 
(A17) 
(A18) 



We generalize the discussion of the parameter f given in SP951 fol- 
lowing the analysis of Giannios et al. (2008). We consider a spher- 
ically expanding shell of cold magnetized jet material of radius R, 
shell thickness A, and Lorentz factor jj, all measured in the lab 
frame. The "spreading radius" of the shell is given by 

^. = 7] A. (Bl) 

The total (isotropic equivalent) energy of the shell is 

E = A7iR^hinmc^y]{\ + a) = M^-^yjC^{\ + cr), (B2) 

where is the rest frame particle number density of the jet mate- 
rial, cr is the magnetization of the material, and M^j is the total rest 
mass of the shell. From E, we obtain the "Sedov length" E and the 
"deceleration radius" Ri^c, 



3E 



1/3 



3Mej 7/1 +tr) 



2/3 



3Mej(l +cr) 



4miim 

1/3 



1/3 



47m I my j 



(B3) 



(B4) 



where iii is the number density of the external ambient medium. 
Substituting for Mq (with R = Rice) in the equation for Rdcc, we 
find that 



i^dec = 3A (1 + (T), 
III 

from which we obtain 

vl/2 



(^dec\ 



— — (1 + 0-) 

Jj "1 



1/2 



(B5) 



(B6) 



Note that ^4 is the number density of the jet ejecta at the moment 
when the shell radius R is equal to R^^.^. 
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